Shear viscosity to entropy density ratio in the Boltzmann-Uehhng-Uhlenbeck model 
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The ratio of shear viscosity [rf) to entropy density (s) for an equilibrated system is investigated in 
intermediate energy heavy ion collisions below 100 A MeV within the framework of the Boltzmann- 
Uehling-Uhlenbeck (BUU) model . After the collision system almost reaches a local equilibration, 
the temperature, pressure and energy density are obtained from the phase space information and 
T]/ s is calculated using the Green-Kubo formulas. The results show that 77/s decreases with incident 
energy and tend towards a smaller value around 0.5, which is not so drastically different from the 
BNL Relativistic Heavy Ion Collider results in the present model. 

PACS numbers: 25.70.-z, 66.20.Cy 



I. INTRODUCTION 

Studying the behavior of nuclear matter under ex- 
treme conditions is one of the most important problems 
in heavy ion collisions. Due to van der Waals nature of 
the nucleon-nuclcon interaction, it is expected that mul- 
tifragmcntation may exhibit features of liquid-gas phase 
transition (LGPT) in intermediate energy heavy-ion col- 
lisions [l|, y]- Evidences of this have been provided 
from various observables, such as the nuclear caloric 
curve, fluctuation, fragment mass distribution and mo- 
ment analysis etc J3l-KL0|. Recent progress on nuclear 
liquid-gas phase transition has been reviewed, es pecially 
for the signals of LGPT in theory and experiment |lll . [l3 |. 

Empirical observation of the temperature or incident 
energy dependence of the shear viscosity to entropy den- 
sity ratio (f?/s) for H2O, He and Ne2 exhibits a minimum 
in the vicinity of the critical point for phase transition 
p^ . Furthermore, a lower bound of rj/s > l/47r ob- 
tained by Kovtun-Son-Starincts (KSS) in certain gauge 
theories is speculated to be valid for several substances 
in nature [14l [l5| . In ultra-relativistic heavy ion collision 
p^l - |2l| . people used the shear viscosity to entropy den- 
sity ratio to study the Quark-Gluon Plasma phase and 
get the minimum value of rj/s, so it is very interesting to 
study shear viscosity or rj/s in intermediate energy heavy 
ion collision [23 - l26| . Unlike the studies on rj/s at rela- 
tivistic energies, there is still very limited investigations 
in intermediate energy heavy-ion collisions. 

In this work, we study the thermodynamic and trans- 
port properties of nuclear reaction and try to see how 
the rj/s evolves with the beam energy or temperature in 
a transport model. We study the equilibration of nuclear 
system within a finite volume using BUU model. In order 
to make the system contain enough number of nuclcons 
in the fixed spherical volume, we choose Au -f- Au system 
in head-on collision (b = Ofm). The system evolves with 



time for long enough time so that it is in freeze-out stage. 
In the final stage of the central collisions, the system 
can be viewed as locally equilibrated. The equilibrium 
in intermediate energy heavy ion collisions can be judged 
by using the temperature and other dynamical variables 
27|. After the system is in equilibrium, we calculate 
the thermodynamic parameters (pressure, energy den- 
sity and entropy density) from phase space information 
of the system. Shear viscosity coefficient is calculated 
from stress tensor fluctuations around the equilibrium 
state using Green-Kubo formula (2^ [23|- Finally, we 
compare rj/s in different incident energy with different 
nuclear equation of state and discuss the results. 

The rest of the paper is organized as follows: In Sec. 
II, we describe the situation of system equilibrium. In 
Sec. Ill, we calculate viscosity coefficient and entropy 
density. Finally, a brief summary and outlook is made in 
Sec. IV. 



II. EQUILIBRATION OF FINITE NUCLEON 
SYSTEM 



We calculate the shear viscosity to entropy density ra- 
tio ij/ s of an equilibrated nuclear system in intermediate 
energy heavy ion collisions using BUU model, which is 
a one-body microscopic transport model based upon the 
Boltzmann equation [23, [30 ■ 



The BUU equation reads [31 1 
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It is solved with the method of Bertsch and Das Gupta 
[33 |. In Eq.([l]), ^^^ and V12 are in-medium nucleon- 
nucleon cross section and relative velocity for the collid- 
ing nucleons, respectively, and U is the mean field potcn- 



tial including the isospin-dcpendcnt term: 



Po Po Po 



(2) 



where po is the normal nuclear matter density; p, p„, 
and pp are the nucleon, neutron and proton densities, 
respectively; r^ equals 1 or -1 for neutrons and protons, 
respectively; The coefficients a, b and a are parameters 
for nuclear equation of state. Two sets of mean field 
parameters are used in this work, namely the soft EOS 
with the compressibility K of 200 MeV (a = -356 MeV, b 
= 303 MeV, a = 7/6), and the hard EOS with K of 380 
MeV (a = -124 MeV, b = 70.5 MeV, a = 2). C^ym is the 
symmetry energy strength due to the density difference of 
neutrons and protons in nuclear medium, here Csym = 32 
MeV is used. 

In this work, we focus on the thermodynamic and 
transport properties of a nuclear system. For this pur- 
pose, we investigate the process of the head-on Au + Au 
collision in a spherical volume with the radius of 5 fm. 
In the following calculation, we show the results with the 
hard EOS except in the case we mentioned. 

First we check the evolution of equilibration situation 
and temperature. The anisotropy ratio, which is a mea- 
sure of the degree of equilibration reached in a heavy-ion 
reaction, is defined as 



Rv 



2 R\\ 



(3) 



where R\\ = {Jp^. + Py) and i?_L = (yT^f) are calculated 

by the momentum of nucleons in the given sphere. As an 
example, the time evolutions of Rp for Au-hAu systems 
within a 5fm-radius sphere at 50 A MeV arc shown in 
Fig. 1(a). When Rp approaches to 1 at around 100 fm/c, 
nuclear system is under equilibrium. 

Time evolution of temperature is also used to judge the 
state of equilibration. Temperature of the system can be 
derived from the momentum fluctuations of particles in 
the center of mass frame of the fragmenting source [Sj] . 
The variance cr^ is obtained from the Q^ distribution 
through 



(Q.') - {Qzf 



(4) 



where Qz is the quadruple moment which is defined by 
Qz = '2-Pz'^ — Px^ ~ Py'^^ ^'^d Px,Py and pz are three com- 
ponents of momentum vector extracted from the phase 
space of BUU model. If the mean equals zero, the second 
term vanishes. Q^ is described by 



{Qz')^ dM^Pz'-p.'-PyTfip). 



(5) 



Assuming a Maxwellian distribution for the momentum 
distribution, i.e. 
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FIG. 1: -Rp (a), temperature (b) and entropy density per nucleon 
(c) as a function of time (after 24 fm/c) for the head-on Au -|- Au 
collision within 5fm-radius sphere at 50 A MeV. 



we can obtain 



{Qz^) = Am^A^T^ 



(7) 



after Gaussian integral, where m is the mass of a nucleon 
and A is the mass number of the fragment. For a nucle- 
onic system, we have A = 1 and can calculate the evolu- 
tion of temperature using this equation. Fig. 1(b) shows 
the temperature's evolution after 25 fm/c, it is seen that 
temperature reaches a maximum around 50 fm/c when 
the system is in the most compressible stage and then it 
starts cool down when the system expands, later on the 
system tends thermodynamic equilibrium. For an equili- 
brated system, the kinetic energy distributions approach 
the Boltzmann distribution as time increases [34|. Af- 
ter the expansion process, the system will approach an 
equilibrate state, so we can investigate the viscosity co- 
efficient and entropy density in system. 

Except temperature, other thermodynamic variables 
can be calculated during heavy-ion collisions. Energy 
density inside a volume with the 5fm radius can be de- 
fined as 



1 

V 



E 

ri<ro 



E^. 



(8) 



where Ei is y/pf + mf, ri is the position of the i-th nu- 
cleon in the center of mass and Tq is the selected radius 
(here we set tq = 5 fm) and pressure can be defined as 



P 



W ^ E,' 

ri<ro 



(9) 



After we get the energy density, pressure and temper- 
ature, entropy density can be calculated by the Gibbs 
formula 



s + P ~ fJ-nP 
T 



(10) 



where /k„ is the nucleon chemical potential and p is nu- 
cleon density of system within the given sphere. In prin- 
cipal, once we have the temperature T and /(p), we can 
fit to a Fermi-Dirac function to extract the chemical po- 
tential. However, we can assume, for simplicity, zero nu- 
cleon chemical potential, or /i„ can be taken around 20 
MeV in the present calculation [35[ . In the following cal- 
culations of entropy and 77/s, we only show the results 
with /x,j = 20 MeV. But we have also checked the results 
with zero chemical potential, this will increase the en- 
tropy about 8% and then lead to the decreasing of -qj s 
about 8%. However, this does not change our conclusions 
of this work. Fig. 1(c) shows the entropy density per nu- 
cleon [s/A) evolves with time after 25 fm/c. It seems the 
entropy density per nucleon reaches a minimum when the 
system is in the most compression stage and rises up in 
the expansion phase, it finally reaches to an asymptotic 
value. From the viewpoint of phase space, the number of 
occupied states is most limited in the high density phase. 
In the high density phase, Fermi blocking forces some of 
the nucleons into higher momentum states but still in a 
spatially confined region, which makes the entropy den- 
sity the minimum. When the system expands, the Fermi 
blocking is reduced, the system cools, but more coordi- 
nate volume is occupied. Therefore it is possible that this 
leads to a non-isentropic expansion until the system is in 
the freeze-out stage. 

Energy density and temperature can be calculated in 
different time when the reaction is going on. Fig. 2(a) 
shows energy density per nucleon versus temperature for 
the studied system after 25 fm/c. Similarly, pressure per 
nucleon increases with temperature is shown in Fig. 2(b). 
From the figure, we can see both energy density and pres- 
sure increase with temperature. This can be understood 
that in a given volume, the increasing of temperature 
reflects stronger thermal motion of nucleons, therefore 
the kinetic energy will give more contribution on energy 
density and pressure. 



III. VISCOSITY COEFFICIENT AND 
ENTROPY DENSITY 

Now let's move on the investigation on transport co- 
efficient of Au-|-Au system in a given volume by the 
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FIG. 2: Energy density per nucleon (a) and pressure per nucleon 
(b) as a function of temperature for the head-on Au -|- Au collision 
within Sfm-radius sphere at 50 A MeV. 



BUU model. Viscosity is one of the transport coefficients 
which characterize the dynamical fluctuation of dissipa- 
tivc fluxes in a medium. Transport coefficients can be 
measured, as in the case of condensed matter applica- 
tions. Also, they should be in principle calculable from 
the first principle. Monte-Carlo simulation for transport 
coefficients is a powerful tool when studying transport 
coefficients using Green-Kubo relations [3^, 123] ■ In high 
energy heavy-ion collisions, calculation of transport co- 
efficients of shear viscosity for a binary mixture |38| , and 
the calculation of coefficient of a hadrons gas has been 
studied [SJ, |33 . The situation of nuclear gas in interme- 
diate energy heavy ion collisions is similar to hadrons gas. 
To study the extended irreversible dynamic processes, we 
use the Kubo fluctuation theory to extract transport co- 
efficients [2^. The formula relates linear transport co- 
efficients to near-equilibrium correlations of dissipative 
fluxes and treats dissipative fluxes as perturbations to 
local thermal equilibrium. The Green-Kubo formula for 
shear viscosity is defined by 



^l = 







(11) 



where T is the equilibrium temperature of the system, t 
is the post-equilibration time (the above formula defines 



t — as the time the system equihbrates and is deter- 
mined by equihbrium time), and {Trij{0,0)iTij{r,t)) is the 
shear component of the energy momentum tensor. The 
expression for the energy momentum tensor is defined by 
the momentum tensor reads 34 



T 



3 "«J ^ » 



Tu{r,t)^ ld'p^f{x,p,t), 



(12) 



where f{x,p,t) is the phase space density of the par- 
ticles in the system. In order to compute an integral, 
we assume that nucleons are uniformly distributed in the 
space. Meanwhile, the isolated spherical volume with the 
radius of 5 fm is fixed, so the viscosity becomes 



V 2 

V= y(7r,j(0) )r^, 
where r^r is calculated by 



(7rij(0)7rjj(t)) ocexp( ). 



(13) 



(14) 



As shown in Fig. 3(a), (7r,j(0, 0)7ry (r, t)) is plotted as a 
function of time for Au -I- Au at 50 A MeV. The correla- 
tion function is damped exponentially with time and can 
be fitted by the Eq. p4)) to extract the inverse slope cor- 
responds as the relaxation time. Fig. 3(b) summarizes 
the relaxation time decrease as the increase of incident 
energy, indicating that the system can approach to equi- 
libration faster at higher incident energy. 

Using the above method, we present the value of rj/s 
as a function of incident energy after the studied system 
has been in equilibrium in Fig. 4. The two sets of nuclear 
equation of state are used. The rj/s value shows a rapid 
fall as the increasing of incident energy up to £■ < 70A 
MeV and then drops slowly to a value close 0.5 when 
E > 70A MeV . Since the BUU equation is one-body 
theory, fragmentation which originates from the fluctu- 
ation and correlation can not be treated in the present 
model. In this case, the phase transition behavior cannot 
be predicted in the BUU model. The continuous drop of 
the ratio of shear viscosity to entropy density does not 
show a minimum at a certain beam energy, which indi- 
cates no obvious phase change or critical behavior in the 
present model. This is a shortcoming of the BUU model 
itself, especially when it is applied to higher beam en- 
ergy. Actually, when we calculate the differential values 
of 77/s versus the beam energy (see the inset of Fig. 4), 
it seems a turning point around E ^ 65A MeV. This 
turning point could indicate the change of dynamical be- 
havior of system, in other words, other mechanisms may 
be needed to be taken into account in the model espe- 
cially at higher beam energies, eg. multifragmcntation 
[40|. Alternatively, lack of minimum 7]/s perhaps shows 
that a behavior with a local minimum of rj/s atphase 
transition temperature might not be universal [41[. In 
the present BUU calculation, all calculated values of r]/s 
are well above the conjectured KSS lower bound of l/47r 
[ij, |l5|. Comparing these values of rj/s for our finite 
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FIG. 3: (a) {nij{0,0)'rTij{r,t)) evolves with time for the head-on 
Au+Au collision in a given 5fm-radius volume at 50 A MeV; (b) 
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Au-f Au collision in a given 5fm-radius volume. 
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FIG. 4: (Color online) rj/s as a function of beam energy for the 
head-on Au + Au collision in a spherical volume with radius of 5 
fm. The inset shows the derivative of rj/s versus beam energy. 



nuclei in BUU model, we see they are not drastically dif- 
ferent either from the RHIC results [13, |42[ or from the 
results of the usual finite nuclei at low temperature from 
the widths of giant vibrational states in nuclei [2g . As 
pointed in Ref. [26[ , it is possible that the strong fluidity 
is a characteristic feature of the strong interaction of the 
many-body nuclear systems in general and not just of the 
state created in the relativistic collisions. Another inter- 
esting point from Fig. 4 is that 77/s shows EOS sensitivity 
in lower beam energy: hard EOS displays larger 77/s than 
the soft one, i.e. larger compressibility of nuclear matter 
can lead to higher 77/s value. 



IV. SUMMARY 

In summary, we studied thermodynamic variables as 
well as viscosity and entropy density for heavy ion col- 
lision after the system tends towards equilibrium in in- 
termediate energy heavy ion collisions in the framework 
of BUU model. The Green-Kubo relation has been ap- 
plied for the nucleonic matter in a central region with a 
moderate volume when the system has been in equilib- 
rium stage for central heavy-ion collisions of Au -|- Au. 
It is found that the ratio of shear viscosity to entropy 
density 77/s decreases very quickly before 70A MeV and 
then drops slowly towards a smaller value of 77/s around 
0.5 at higher beam energy in Au-f Au system. The 77/s 
are not drastically different either from the RHIC results 
or from the results of the usual finite nuclei at low tem- 
perature. However, no obvious minimum 77/s value oc- 



curs within the investigated energy range. This may re- 
fiect that no liquid-gas phase transition behavior is dis- 
played in the present model due to the shortcoming of 
the model itself which lacks dynamical fluctuation and 
correlation. Relating the shortcoming, the equilibrium 
temperature could be a little higher than other model 
(QMD, SMM) which considers fragment formation [i^ 
and while the shear viscosity and entropy density might 
be influenced by the cluster formation. Therefore, other 
models which can incorporate liquid gas phase transition 
should be checked for the shear viscosity and entropy 
density. For instance, it will be very interesting to use 
quantum molecular dynamics-type model to check if a 
minimum of 7;/s will occur around liquid-gas phase tran- 
sition. The work along this direction is in progress. Of 
course, experimental studies on shear viscosity is more 
important to demonstrate the relation of ?7/s and liquid- 
gas transition point. 
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